General relativistic statistical mechanics 



(N 

o 

(N 

> 
O 

OV 



> 
IT) 

o 
q 

ON' 

o: 



X 



Carlo Rovelli 

Centre de Physique Theorique de Luminy, Aix-Marseille University, F-13288 Marseille, EU 

Understanding thermodynamics and statistical mechanics in the full general relativistic context is an 
open problem. I give tentative definitions of equilibrium state, mean values, mean geometry, entropy 
and temperature, which reduce to the conventional ones in the non-relativistic limit, but remain valid 
for a general covariant theory. The formalism extends to quantum theory. The construction builds 
on the idea of thermal time, on a notion of locality for this time, and on the distinction between 
global and local temperature. The last is the temperature measured by a local thermometer, and is 
given by kT = hdr/ds, with k the Boltzmann constant, h the Planck constant, ds proper time and 
dr the equilibrium thermal time. 



I. INTRODUCTION 

Thermodynamics and statistical mechanics are power- 
ful and vastly general tools. But their usual formulation 
works only in the non-general-relativistic limit. Can they 
be extended to fully general relativistic systems? 

The problem can be posed in physical terms: we do not 
know the position of each molecule of a gas, or the value 
of the electromagnetic field at each point in a hot cavity, 
as these fluctuate thermally, but we can give a statistical 
description of their properties. For the same reason, we 
do not know the exact value of the gravitational field, 
which is to say the exact form of the spacetime geome- 
try around us, since nothing forbids it from fluctuating 
like any other field to which it is coupled. Is there a 
theoretical tool for describing these fluctuations? 

The problem should not be confused with thermody- 
namics and statistical mechanics on curved spacetime. 
The difference is the same as the distinction between 
the dynamics of matter on a given curved geometry ver- 
sus the dynamics of geometry itself, or the dynamics of 
charged particles versus dynamics of the electromagnetic 
field. Thermodynamics on curved spacetime is well un- 
derstood (see the classic [l|) and statistical mechanics on 
curved spacetimes is an interesting domain (for a recent 
intriguing perspective see Q). The problem is also dis- 
tinct from "stochastic gravity" [!, H| , where metric fluc- 
tuations are generated by a Einstein-Langevin equation 
and related to semiclassical effects of quantum theory. 
Here, instead, the problem is the just the thermal behav- 
ior of conventional gravity^] 

A number of puzzling relations between gravity and 
thermodynamics (or gravity, thermodynamics and quan- 
tum theory) have been extensively discussed in the litera- 
ture [5l4l4| . Among the most intriguing are probably Ja- 
cobson's celebrated derivation of the Einstein equations 



1 One may ask whether equilibrium can ever be reached, given the 
gravitational instabilities and long thermalization times. The 
question is legitimate, but doesn't authorize us evading the issue 
of what equilibrium means: for the question itself to even make 
sense, and because we are always concerned only with approxi- 
mate equilibrium in nature, gravity or not. 



from the entropy- area relation |15l[T^ |. and Penrose Weil- 
curvature hypothesis [l?], HI]- These are very suggestive, 
but perhaps their significance cannot be evaluated until 
we better understand standard general covariant thermo- 
dynamics. 

One avenue for addressing the problem is perturbation 
theory. Another is restricting to asymptotic flatness and 
observablcs at infinity |19l - t2l| . Although useful in spe- 
cific contexts, these roads are incomplete, because they 
miss the core issue: understanding if temperature has a 
meaning in the bulk of spacetime in a strong field regime. 
What do we mean when we say that near a cosmological 
singularity temperature is high? For the moment we do 
not have a definition of temperature that makes sense 
where the metric might fluctuate widely. 

A step towardsgeneral covariant statistical mechan- 
ics was taken in [22 . [H| and extended to quantum field 
theory in 24]. The notion introduced in these papers 
is thermal time. This is meant to address the basic dif- 
ficulty of general relativistic statistical mechanics: in a 
generally covariant theory, dynamics is given relationally 
rather than in terms of evolution in physical tim^l, conse- 
quently the canonical hamiltonian vanishes, and without 
a hamiltonian H it is difficult to even start doing statisti- 
cal physics. The idea of thermal time is to reinterpret the 
relation between Gibbs states (p oc e~° H ) and time flow 
(generated by H): instead of viewing the Gibbs states 
as determined by the time flow, observe that any generic 
state generates its own time flow. The time with respect 
to which a covariant state is in equilibrium can therefore 
be read out from the state itself. The root of the tempo- 
ral structure is thus coded in the non commutativity of 
the Poisson or quantum algebra [24], HH . 

Since any state is stationary with respect to its own 
flow, the problem left open is characterizing the states 
that are in physical equilibrium. Here we consider a so- 
lution: equilibrium states are those whose thermal time is 
a flow in spacetime These, we suggest, are the proper 
generalization of Gibbs states to the general covariant 
context. 



2 For a discussion of this crucial point see the Appendix and Chap- 
ter 3 of l2. r | . in particular Section 3.2.4. 



This problem is considered also in [27142 
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This step allows temperature to be defined, following 
the intuition in [53, H(| : the temperature measured by a 
local clock is the ratio between thermal time and proper 
time. This yields immediately the Tolman-Ehrenfest law 
[3ll [32[, which correctly governs equilibrium tempera- 
ture in gravity. Entropy and free energy can be defined 
and we obtain the full basis of generally covariant ther- 
modynamics. The construction extends to the quantum 
theory. 

The result is a tentative set of equations that generalize 
conventional thermodynamics and statistical mechanics 
to classical and quantum general-covariant systems. 

We use units where the Boltzmann constant k and the 
Planck constant ft are set to unity. We have tried to keep 
the main text brief, confining background material to a 
detailed Appendix. The reader is urged to start from the 
Appendix unless the language and the background ideas 
of the text are already familiar. Equations in the paper 
are to be understood locally in phase space, namely on a 
chart where suitable regularity conditions are satisfied to 
avoid singular or degenerate behavior. A finer analysis 
will make sense after the basic conceptual structure is 
clear. 



II. GENERAL COVARIANT GIBBS STATES 
A. Thermal time 

Let £ be a symplectic space, whose physical interpre- 
tation is the extended phase space of a general covariant 
theory (see the Appendix for notation and details.) Let 
C be a submanifold of £ , representing the surface where 
the constraints of the theory (which code the full dynam- 
ics) are satisfied. The symplectic form a of £ induces a 
presymplectic structure on C, whose null directions can 
be integrated to define the gauge orbits o. The space T of 
these gauge orbits, which is the physical phase space of 
the theory, is again a symplectic space, with symplectic 
form ui. It is in 1-to-l correspondence with the space of 
the solutions of the field equations, modulo gauges. A 
statistical state p is a positive function on T normalized 
with the respect to the Liouville measure 

The hamiltonian vector field X defined by 

p u(X) = dp (2) 

generates a flow a T in T called the thermal flow; its gen- 
erator 

ft = - In/j (3) 



is called the thermal hamiltonian and the flow parameter 
t is called thermal tim^3 (22j. 



B. Local thermal time 

Consider a general covariant theory that includes gen- 
eral relativitjQ and assume physical 3d space S to be 
compact, with the S3 topology. The space £ can be 
coordinatized by the 3d Riemann metric tensor q of E, 
the matter fields tp, and their respective conjugate mo- 
menta (p, 7r); these quantities are fields on E, namely 
functions from E to a target space {q,tp,p,ir) ■ E — > V. 
An orbit o determines a solution of the field equations 
and therefore in particular a pseudo-riemannian mani- 
fold (M, g) a , A point in o determines a spacelike Cauchy 
surface cf> : S3 — > (M,g) a , having the given induced met- 
ric q and extrinsic curvature p. In particular, a foliation 
4> T : S3 — > (M,g) ,T e R of (M,g) a corresponds to a line 
on the orbit. 

Consider now a real function T on V. This deter- 
mines a local function (which we indicate with the same 
letter) on £ , namely a map T : £ x E — > R given by 
f((q,p,<p,*l,x) = r(g(x),p(x) l¥ >(x),7r(x)),x G E. The 
coordinate T(x) on £ plays a role of "multi-fingered time" 
in what follows. If the equation 

f (x) = t, t G R (4) 

defines a foliation of (M, g) Q (on a given region of phase- 
space) we say that T(x) is a "local time". The pa- 
rameter of the foliation defines then a time coordinate 
t : [M, g) ~ > R on spacetime. The simplest example 
is if the matter fields include a scalar field T that grows 
monotonically in spacetime (for the given region of phase 
space): then the value of the field defines a time coordi- 
nate. 

If there are canonical coordinates T and Q l on £ , with 
respective momenta Pf and Pi, such that 

P f (x) = -ft(Q I (x),P l (x)) (5) 

on C, then T(x) defines a deparametrization of the theory 
in the following sense: the hamiltonian 

h= [ d 3 x ft(Q*(x),P(x)) (6) 



4 So defined, r has the dimensions of an action, as it is conjugate 
to a dimensionless quantity. It can be made dimensionless by 
multiplying the r.h.s.of and J3J by h. This is a bit artificial 
in the classical theory, but will be natural in the quantum theory. 

J We systematically disregard at this stage the difficulty of defining 
the Liouville measure that defines the integral JT} in the case 
of field a theory. This is because the issue should properly be 
addressed in the quantum context, where I will be a bit more 
precise. 
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evolves geometry and matter fields along the foliation <p T . 
Notice that h is constant along the orbits it generates. 
We can therefore associate its value to each o and obtain 
in this manner a function h onT. (Weaker cases are also 
of interest, in particular the case relevant in cosmology 
where 

P f (x) = -/(f(x)) /i(Q'(x),P,(x)); (7) 

which describe a system with temperature varying in 
time: see [23|.) 

Let us now come to the first main notion that we in- 
troduce in this paper. We say that a statistical state p 
on r is a "Gibbs state" if there is a local time T{x) with 
a local hamiltonian h of the form ([6]) (or ([7]l) satisfying 
([3]) up to an additive constant^]. 

If this is the case, the thermal time r generated by p is 
precisely the foliation time r, and therefore thermal time 
has a geometrical interpretation as a flow in spacetime. 

C. Nonrelativistic limit 

The definition above is a generalization of the conven- 
tional definition of Gibbs states. To see this, recall that 
for a Hamiltonian system with phase space To, canon- 
ical coordinates (<?,p) and Hamiltonian H = H{q 1 p), a 
Gibbs state is a state of the form pp = Z^ 1 {/3) e~P H with 
Z((3) = J To e~^ H . The general covariant formulation of 
this system is defined on the extended phase space £ 
with canonical coordinates (t,p t ,q,p) and the constraint 
C = pt + H(q,p). The constraint surface is coordina- 
tized by (t,q,p) and the orbits are given by (t,q(t),p(t)) 
where q(t) and p(t) are the solutions of the Hamilton 
equations. The space T of these orbits is isomorphic to 
To (but not canonically so, until a t = to is chosen) via 

(q,p) = q(to),p(to)- 

A time function on £ is provided by t = t/ ft, whose 
conjugate momentum is p T = f3p t , which satisfies the 
requirement that the constraint can be expressed in the 
form (O, namely p T = —h(q,p) 1 where h = j3H. The 
hamiltonian, being constant on each orbit, is well defined 
on r, therefore pp is a function on L, namely it is a 
statistical state in the covariant sense. It is immediate 
to see that it satisfies (J3J. In other words, the Gibbs 
state picks out the coordinate t from £ , where this was 
confounded with the other variables. 

Observe now that the temperature T = j is equal to 
the ratio 

r = I (8) 

between the thermal time r, namely the parameter of the 
evolution generated by the logarithm of the Gibbs state, 
and the physical time t. This characterization of temper- 
ature can be extended to general covariant systems. 



D. Mean values, mean geometry 
and local temperature 

Consider a family A of functions A on T. Let the mean 
value of A on the state p be 

A = J r A P- ( 9 ) 

The thermal time flow a T acts on these functions by 
A(t)(s) = a T (A)(s) = A(a- T (s)), s e T, which satis- 
fies dA/dr = {A.h}. Since p is clearly invariant under 
the flow, so are the mean values, but 

f AB {r) = Ja(t)B p. (10) 

is in general a non trivial function and describes temporal 
correlations in the state. Define the mean geometry g (if 
it exists) of a state p for an observable family A as a 
spacetime (M, g) with a foliation <fi T such that 

A{r)=A{cj>-\g)). (11) 

Since A(t) is r independent, it follows that (Af, g) is 
stationary under the flow defined by <f> T . Therefore £ = 
|: is a timclikc Killing field on (M,g). The norm of 
£ is ds/dr namely the ratio between the local flow of 
proper time and thermal time. The equivalence principle 
therefore compels us to define the local temperature by 
the local version of ©, namely 

T(x) = \S(x)\- 1 , xeM (12) 

from which the Tolman-Ehrenfest law [3ll [32j 

T(x) |£(x) | = constant (13) 

that governs the spacetime variation of temperature at 
equilibrium in gravity, follows immediately^ In station- 
ary coordinates (r, x) , the temperature is the inverse of 
the Lapse function, since ds 2 — N 2 dr 2 . 



E. Partition function 
and global temperature 

If p is a Gibbs state, we can obtain another Gibbs state 
by exponentiating it with a constant (3 and multiplying it 
by a P dependent factor that preserves the normalization: 
pp = Z _1 (/3) pP . The effect of this exponentiation is to 
scale the thermal time globally, and therefore to scale the 
temperature globally. Therefore the global temperature 
is defined with respect to a reference Gibbs state. Having 



The constant has no effect on the dynamics and we set it to zero 
by redefining Pf. 



7 A suggestion in this direction was in [27J . The intriguing relation 
between (JSJl and the Tolman law was pointed out in 13011 . 
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a one-parameter family of Gibbs states allows us to define 
the partition function 

Z(J3) = Jff t . (14) 

The entropy of the state can be obtained as usual from 

S{p) = -Jpphipp (15) 

and from this we can derive in a few steps the standard 
thermodynamical relation 

S = - J p, 3 {fi\n P -\nZ) 

= PE + lnZ (16) 

where E is the mean value of the energy h = — In p. The 
global temperature fj of the state should not be con- 
fused with the local temperature T(x), which is space- 
dependent. Also, the local temperature is defined di- 
rectly by a single statistical state (if a mean geometry 
exists), while the global temperature is only defined rel- 
ative to another Gibbs state taken as reference. 

In the following section, we extend this structure to 
the quantum theory. 

F. Quantum theory 

Let K. be the unconstrained Hilbert space of a general 
covariant theory and % its physical Hilbert space (the 
"space of solutions of the Wheeler-deWitt equation"). 
General covariant quantum mechanics is well defined by 
these structures (See the Appendix, and, more in detail, 
Section 5.2 of |||.) 

A quantum statistical state is a trace-class operator p 
on H such that trp =1. Its entropy is S — — trfplnp]. 
Let A be an observable algebra formed by self-adjoint 
operators A on T~L. Then p defines a state on this algebra 
by 

p(A) = to[Ap] (17) 

and3 the Tomita theorem provides a flow a T : A — > A on 
the observable algebra. This is the thermal-time flow in 
the quantum theory 24] . If there is a local hamiltonian h 
and a (now dimensionless) conjugate "time" observable 
r that in the classical theory reduces to the quantities 
defined in the previous section, and generates an evolu- 
tion 

a T (A) = ei hT Ae-% hr , (18) 



Taking A to be a von Neumann algebra, namely a * -algebra 
of bounded operators closed in the weak operator topology and 
including the identity. 



then we say that p is a Gibbs stateH The Tomita flow of 
p satisfies the KMS condition (see, for instance, [33]) 

fAB(r) = fBA(-T + 2TTi) (19) 

for any two observables A and B, where 

Jab(t) =p(a T (A)B). (20) 

A thermal state p p = Z~ x {fi) p 0/27r satisfies the KMS 
condition 

fAB(r) = f BA (-t + ip) (21) 

with respect to the flow generated by p. 

The notion of mean geometry can be extended to the 
quantum theor-J^l by defining (M,g,(j) T ) (if it exists) as 
the mean geometry of the state p with respect to a given 
observable algebra A if 

A(r)^p(a T (A)))=A(^ 1 (g)). (22) 

The local temperature T(x) is defined by the norm of 
the killing field of the mean geometry, and is therefore 
a semiclassical concept. Restoring physical units, local 
temperature is given on the mean geometry by 

„, , h dr , . 

T ™ = k ™ 

where h is the Planck constant and k is the Boltzmann 
constant. 

Notice that (|23|) gives the Unruh temperature [Icj of a 
quantum field theory on Minkowski space, if ds is the 
proper time along the accelerated observer trajectory 
and t is the dimensionless parameter of the Bisognano- 
Wichman flow U(t) — e lrK / 2 ^ i where K is the boost 
generator, which is the Tomita flow of the vacuum state 
on the Rindler- wedge observables [H, EH ■ 

This suggests that the Unruh effect should affect the lo- 
cal temperature of an observer accelerated on a mean ge- 
ometry, also in the context of the full generally-covariant 
statistical mechanics of the gravitational field. If a mean 
geometry has a Killing horizon, where the norm of £ be- 
comes singular, then the local temperature (|23j) diverges 
on the horizon. The divergence of the temperature is a 
high-energy, namely a short-distance phenomenon, there- 
fore we can consider it in a region of spacetime small 
with respect to the local curvature of the mean geome- 
try, namely as a locally flat-space phenomenon. As such, 
it must be determined by the Unruh temperature. An 
explicit example of a statistical state where this happens 
has been discussed in [42|, |43[ • An Unruh temperature in 
the vicinity of the horizon of a black hole is red-shifted 
by the Tolman relation (fT3|) precisely to Hawking's black 
hole temperature at infinity. 



Since space is compact, the usual difficulty of hamiltonian quan- 
tum field theory with thermal states which historically gave rise 
to algebraic quantum field theory, is not there, since energy does 
not diverge on thermal states. 

The idea of mean geometry is implicit in contexts where covariant 
quantum states of gravity are associated to a classical geometry 
34- 39]. 
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III. CONCLUSION 



APPENDIX 



We have extended the machinery of statistical ther- 
modynamics to the general covariant context. The new 
concepts with respect to conventional statistical mechan- 



A. Classical theory 

1. Mechanics 



1. The statistical state is defined on the space of the 
solution of the field equation. 

2. Each statistical state defines a preferred time flow, 
called thermal time. 

3. A statistical state whose thermal time flow has a 
geometrical interpretation, in the sense that it can 
be reinterpreted as evolution with respect to a lo- 
cal internal time, defines a generalized Gibbs state, 
with properties similar to the conventional equilib- 
rium states. 

4. For such states, it is possible to define the relative 
global temperature between two states. 

5. A mean geometry is a stationary classical geome- 
try with a timelike killing field and a time foliation, 
such that the value of a suitable family of observ- 
ables reproduces the statistical expectation values 
of these observables in the statistical ensemble. 

6. If a mean geometry exists, a local temperature is 
defined. Local temperature is the ratio between 
proper time and thermal time on the mean geome- 
try: 



T(x) 



H dr 
k ds ' 



(24) 



It yields immediately the Tolman law. 



This construction reduces to conventional thermodynam- 
ics for conventional Hamiltonian systems rewritten in a 
parametrized language. 

Examples, extension of the formalism to the bound- 
ary formalism (4^ - [46j j . which is the natural language for 
quantum field theory in the generally covariant context, 
and applications to horizon thermodynamics, and in par- 
ticular to the local framework defined in [47| and the 
derivation of black hole entropy in loop quantum gravity 
in 42[, will be considered elsewhere. 



I thank Alejandro Perez for the crucial suggestion to fo- 
cus on the locality of the hamiltonian, Ed Wilson- Ewing 
for pointing out the relevance for cosmology of the weaker 
notion of equilibrium, captured by (JTJ), Eugenio Bianchi 
for numerous discussions on this subject, Simone Speziale 
and Pierre Martinetti for several helpful comments. 



A conventional hamiltonian system is defined by a 2N 
dimensional phase space Tq and a hamiltonian H. The 
phase space is a symplectic space, namely a manifold 
equipped by a non-singular closed symplectic two form 
uj. Locally, we can always choose coordinates (q l ,Pi) on 
T such that 



u) = dq 1 A dpi 



(25) 



(summation understood). Having a symplectic two form 
is the same as having Poisson brackets. H is a scalar 
function on V . Every function / on a symplectic space 
defines a vector field Xj on the space, defined by 



,(X f ) = -df, 



(26) 



where the l.h.s is the action of a differential two-form on 
a vector, which gives a one form, and the r.h.s is the 
differential of /. In turn, a vector field defines a flow 
a t : Tq — > Tq, t £ R, namely a continuous one-parameter 
group of automorphisms of Tq into itself, related to X by 



da t 
~~dt 



X 



(27) 



t=o 



The Poisson bracket between two functions A and B on 
r is defined by 



{A,B}=X B {A) = -X A {B). 



(28) 



The flow of the hamiltonian is the time flow, namely the 
evolution in time of each point of Tq. Explicitly, the 
hamiltonian vector field of H is easily seen to be 



X = 



dH d dH d 

dpi dq 1 dq 1 dpi ' 



(29) 



so that the time flow is determined by the Hamilton equa- 
tions 



dq l {t) 
dt 



dH 

dp i ' 



dpj(t) 
dt 



dH 



(30) 



which show that this geometric construction is equivalent 
to hamiltonian mechanics. An observable A is a real 
function on Tq. The time evolution of an observable is 
defined by A(t) = Ao a t and satisfies 



dA(t) 
dt 



= {A,H}. 



(31) 



Let F be the space of the solutions of the equation 
of motion (q 1 (t) , Pi(t)) . This is a finite dimensional space 
with is isomorphic to r , but not canonically isomorphic. 
A specific isomorphism is obtained by choosing a value 
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for the time parameter t. Then the isomorphism between 
T and T is given by (q l ,Pi) — {q l (t ) , Pi(t j) ■ Thanks to 
this isomorphism, T has a symplectic structure as well 
(independent from to). 

For instance, the solutions of the dynamics of a 
harmonic oscillator have the form (q(t) = Asm(u)t + 
<f>), p(t) = mujAcos(ujt + 4>)). Therefore T is co- 
ordinatized by A and <j). A map between T and T 
is obtaining choosing for instance t = 0, which gives 
(q = Asincj), p — mu; A cos <f>) and therefore the sym- 
plectic form on T is 

a = -mu AdAAdcj). (32) 

An equivalent formulation of the dynamics, called the 
presymplectic formulation, can be given on the 2N + 1 
space C = To x R, with local coordinates (q l ,pi,t) 
equipped with the two-form 

J = dq l Adp l -dH{q l ,p i ) Adt. (33) 

This two- form is degenerate, namely has a null direction 
(since the space is odd dimensional). That is, there exist 
a vector field X', determined up to scaling, such that 

oj'(X') = 0. (34) 

It is immediate to sec that this vector field is proportional 
to 

X ' = lT t +X <»> 

and its integral lines (called the orbits of u') are precisely 
(the graphs of the) physical motions (t,q l (t),pi(t)). Let 
r be the space of these orbits and tt the projection that 
sends each point of C to the orbit to which it belongs. 
T carries a symplectic two-form a uniquely characterized 
by the fact that its pull back to C by tt is u/. (A pull 
back is degenerate in the directions of the orbits.) The 
symplectic space (r,cr) is clearly the same as the one 
constructed above. 

The equivalence between the conventional hamilto- 
nian and the presymplectic formulation, is almost com- 
plete. The reason for the "almost" is subtle, interesting, 
and at the core of the problem discussed in this paper. 
Given a hamiltonian system (Tq,uj, H), we can immedi- 
ately construct its corresponding presymplectic formula- 
tion (C,u'). But the opposite is not true, since we need 
to know which one of the variables on C is the time vari- 
able, in order to do so. In other words, the presymplectic 
formulation leads to the same relations between the vari- 
ables (t,q l ,pi) as the hamiltonian one, but without spec- 
ifying which of these variables is to be recognized as the 
time variable. The difference is the same as the differ- 
ence between giving a function y(x) or its parametrized 
form (y(s), x(s)): in the first case x is singled out as the 
independent variable, in the second case it is not. 



2. General covariant mechanics 

Systems like general relativity, or a single free relativis- 
tic particle, are defined in the covariant language by a 
Lagrangian that leads to a vanishing canonical Hamilto- 
nian. Equivalently, they are defined by equations of mo- 
tion that are gauge invariant under a re-parametrization 
of the evolution coordinate. The Legendre transform of 
the Lagrangian of these systems defines a phase space 
with constraints, and the dynamics is coded in the con- 
straints. Let £ denote this phase space (to distinguish 
it from the phase space of a conventional system, since 
it has a different physical interpretation) and let C de- 
note the subspace of £ where the constraints vanish. £ 
is a symplectic space with symplectic form lj. Its restric- 
tion to C is a presymplectic two-form u/ (the pull back 
of bj under the embedding i of C in £), which is degen- 
erate in the directions of the hamiltonian vector fields of 
the constraints themselves. The space of the orbits T is 
again a symplectic space carrying a symplectic two-form 
cr, uniquely characterized by 

ixLU = Uj' = 7T*<7. (36) 

where 

£ *±- C T. (37) 

The presymplectic constraint surface (C,a/) defines the 
dynamics precisely as in the presymplectic formulation of 
the hamiltonian dynamics described above. Notice that 
it defines all the physical correlations among dynamical 
variables, without specifying one of these as the indepen- 
dent time variable. The distinctive feature of the general 
covariant systems is therefore to define dynamics as a 
"democratic" correlation between variables instead of as 
evolution with respect to a singled out independent vari- 
able. 

A simple example is provided by the dynamics of a 
free relativistic particle. The extended phase space £ is 
8-dimensional, with coordinates (x^,^) and co — dx^ A 
dp^). The constraint surface C is given by p 2 = to 2 . The 
orbits are given by 

^(t) = ^t + <. (38) 

TO 

And there is a six dimensional space of these. Each or- 
bit determines a correlation between observables. For 
instance it determines the relation between different co- 
ordinates on Minkowski space. Notice that all this is 
Lorentz invariant. Notice also that this canonical formu- 
lation never specifies one particular Lorentz time as the 
preferred one. To obtain a conventional Hamiltonian for- 
mulation we have instead to select a Lorentz frame and 
choose one variable, say x° (as opposed to 2:° = A°x M 
where A is Lorentz matrix) as the time variable. Then 
this determines a Hamiltonian H = ^Jp 1 + to 2 , which 
generates the same motions, but in a non-manifestly 
Lorentz invariant language. 
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Notice that any Gibbs state 



(39) 



breaks Lorentz invariance and selects a preferred Lorentz 
time. Physically, this is the specific Lorentz-time flow 
with respect to which a given gas of relativistic particles 
is in equilibrium. 

Notice that it is somewhat misleading to state that the 
full dynamics of a generally covariant system is entirely 
captured by the physical phase space T and all functions 
on r, because this would be like saying that the dynamics 
of a harmonic oscillator is captured by writing down the 
phase space coordinated by A and 0, and all functions 
of A and <fr. If we do so, we loose track of the fact that 
the harmonic oscillator is characterized by the oscillat- 
ing variable q(t)\ The dynamics of a generally covariant 
system is not just described by T and the family of all 
functions on V. We also need to give explicitly the em- 
bedding of each orbit in C or, equivalently, in £. In the 
case of the relativistic particle, for instance, the dynamics 
is not just the specification that the physical space is six 
dimensional: it is also the information that each point of 
this space determines a timelike line in Minkowski space, 
namely a correlation between quantities on C. In this 
context, such quantities are called "partial observables" 

H. 



3. Statistical mechanics 

The symplectic form defines a volume- form on To, ob- 
tained by taking N times the wedge product of u> with 
itself. This defines an integral on To, which we indicate 
simply without measure notation. A statistical state is a 
real non-negative function p on Tq normalized as 

p=l. (40) 

Its entropy is defined by the Shannon expression 

S = -Jp\np. (41) 

The mean value of an observable A in the state p is de- 
fined by 



.4 



Ap. 



(42) 



The mean value of A{t) can be equally obtained as the 
mean value of A on the state p(t) which satisfies 



dp(t) 
dt 



{P,H}. 



(43) 



An equilibrium Gibbs state, is a particular statistical 
state of the form 



where j3 = 1/kT is a positive real number and T is the 
temperature. It is immediately clear that a Gibbs state 
is time independent and the mean value of all observables 
in a Gibbs state arc time independent. Nontrivial time 
correlations can nevertheless be defined from quantities 
like 



f AB (t) = / A(t)B(0)p. 



(45) 



The proportionality factor in (|44p is determined by the 
normalization condition: 



1 



P 



-PH 



where 



Z{fi) = I e~ m = e-? F 



(46) 



(47) 



is called the partition function, and F is called the free 
energy. It follow immediately from the definitions and a 
short calculation that the mean value E of the energy is 
given by 



and 



lrflnZ 



S = j3E- f3F. 



(48) 



(49) 



these are the basic thermodynamical relations for the 
Gibbs states. 



4- General Covariant Statistical Mechanics 

Here a condense the results of this paper. A statistical 
state is a normalized positive function on the physical 
state space. It determines a thermal flow with generator 
X defined by 



p oj{X) = dp. 



(50) 



p oc e 



-PH 



(44) 



the generator of this flow is the (state dependent) thermal 
hamiltonian h — — In p and the thermal time r is the 
parameter of this flow. For a conventional Gibbs state 
in a non-generally-covariant system, temperature is the 
ratio between thermal time and geometrical time. 

In a gravitational field theory, if h is local, then it 
defines a flow in spacetime, and a preferred foliation of 
the mean geometry. The local temperature, which satis- 
fies the Tolmann relation, is the local ratio between the 
spacetime flow and proper time. 



B. Quantum theory 

1. Quantum Mechanics 

A conventional quantum system is defined by a Hilbert 
space T-L and a family A of observables A, self-adjoint 
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operators on "H, which in particular includes a Hamilto- 
nian H. The Hamiltonian generates a unitary flow on % 
by the one-parameter group of unitary transformations 
U (t) = er %Ht in the Schrodinger picture and a flow on 
the observables by A(t) = U(—t)AU(t) in the Heisen- 
berg picture. In the Schrodinger picture, the Hilbert 
space Ti corresponds to the phase space Tq at a given 
time; while in the Heisenberg picture the Hilbert space 
H corresponds to the phase space T of the solutions of 
the equations of motion. The expectation value of an 
observable in the state tp € % is given by A — (i/jAip) , or 
cquivalcntly by 

A = tv[Ap] (51) 

where 

p=|V><# (52) 

The eigenvalues of A determine the quantization, namely 
the possible outcomes of a measurement, of A and tran- 
sition probabilities between such measurement outcomes 
are determined by the matrix elements of U(t) in the 
observable's eigenbasis. 



equation" , and a family of observables A, B on JC called 
"partial observables" . 

The eigenvalues of the partial observables determine 
the quantization, namely the possible outcomes of a mea- 
surement [491 - 45 ll ] . and transition probabilities between 
such measurements' outcomes are determined by the ma- 
trix elements 

(q\PW) (56) 

of the (generalized) projection 

P:JC->H, (57) 

in the observables' eigenbases \q) in K (see (25|, Chapter 
5 and (52J). A specific example of a definition of these 
transition amplitudes, finite to all orders, is provided by 
covariant loop quantum gravity [53]. The quantum me- 
chanics of generally covariant systems can therefore be 
well defined without the need of specifying a time vari- 
able. 



4- General Covariant Statistical Quantum Mechanics 



2. Quantum Statistical Mechanics 

A statistical state p is a trace-class operator on H nor- 
malized by 

tr[p] = 1. (53) 

The mean value of an observable in such a state is still 
given by (|5Tj) . The states of the form (f5"2"]) satisfy p 2 = p, 
are called "pure" , and their conventional physical inter- 
pretation is that the probabilistic nature of the uncer- 
tainty in the predictions derived from them is not due 
to our ignorance, but to irreducible intrinsic quantum 
uncertainty. The von Neumann entropy of the state p 

5 = -tr[plnp] (54) 

vanishes on pure states. A Gibbs state is a state of the 
form p oc e~° . The partition function is the inverse of 
its normalization, namely 

Z(fi)=\x[-eP H }. (55) 

Again the basic thermodynamical relation (|4T))) follows in 
a few steps from these definitions. 

3. General Covariant Quantum Mechanics 

A generally covariant quantum system is defined by an 
extended Hilbert space fC, a (possibly generalized^!) sub- 
space H, the "space of solutions of the Wheeler-deWitt 



Namely a subspace of 5 in a Gelfand triple S D K, D S. 



The thermal-time flow of a generally covariant statis- 
tical quantum state p is defined by its Tomita flow. This 
can be constructed as follows. The expectation value of a 
statistical state p on the algebra A of the gauge-invariant 
observables a define a state on this algebra. Assuming 
A to be a C*-algebra, the GNS construction defines a 
Hilbert space T-L where observables are represented by op- 
erators and ip is a vector (even if -0 is a statistical state) . 
Let then S be the operator defined by Sayb — a*ip. It is 
always possible to write S in the form S = Je h / 2 , where 
J is antinunitary and e h is self-adjoint. The Tomita flow 
on the algebra is then defined by 

a t a = e- lth a e tth (58) 

and the Tomita theorem states that this is a one- 
parameter group of automorphisms of the algebra. 

To understand what is going one, start from a nor- 
mal quantum field theory. Pure states are vectors in 
Fock space. Mixed states are density matrices, namely 
trace class operators p on Fock space. These form an 
Hilbert space, which we can call T-L: notice that a sta- 
tistical state p is now represented by a vector in this 
Hilbert space, for which a convenient notation is \p). If 
a is an observable on Fock space, we can represent it on 
% as a\p) — \ap), which is again trace class. If p is a 
Gibbs state for a Hamitonian H at inverse temperature 
j3, namely p — e~^ H , then a straightforward calculation 
shows that J\k) = \k*) and e h \k) = \ e -^ H ke~ pH ) sat- 
isfy the definition of S. Therefore the Tomita flow of 
the Gibbs state is precisely the time flow scaled by the 
temperature: ata = e~ lt ^ H ^ ae %tly ^ H ^ . In other words, 
the Tomita relation between a state and a flow is the 
quantum field theoretical version of the classical relation 
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between a state on phase space and its Hamiltonian flow. 
The operator J flips creation and annihilation operators 
of the quanta over the thermal state, and therefore codes 
the split between positive and negative frequencies. For 
a more detailed discussion, see [24J. Time flow is fully 
coded into the statistical state. The local relation be- 
tween thermal time dr proper time dt and temperature 



T is given by equation (|24l) . 

Thanks to Hal Haggard for a careful reading of the 
manuscript and useful comments. 
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